Effect of short- and long-range scattering in the conductivity of graphene: Boltzmann 

approach vs tight-binding calculations 
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We present a comparative study of the density dependence of the conductivity of graphene sheets 
calculated in the tight-binding (TB) Landauer approach and on the basis of the Boltzmann theory. 
The TB calculations are found to give the same density dependence of the conductivity, cr ~ n, 
for short-range and long-range Gaussian scatterers. In the case of short-range scattering the TB 
calculations are in agreement with the predictions of the Boltzmann theory going beyond the Born 
approximation, but in qualitative and quantitative disagreement with the standard Boltzmann ap- 
proach within the Born approximation, predicting a — const. Even for the long-range Gaussian 
potential in a parameter range corresponding to realistic systems the standard Boltzmann predic- 
tions are in quantitative and qualitative disagreement with the TB results. This questions the 
applicability of the standard Boltzmann approach within the Born approximation, commonly used 
for the interpretation of the results of experimental studies of the transport in graphene. 

PACS numbers: 73.63.-b , 72.10.-d , 73.63.Nm. 73.23.Ad 



I. INTRODUCTION 

Understanding factors that affect conductivity of 
graphene represents a fundamental task of great impor- 
tance in view of possible application of graphene-based 
devices for electronics and optoelectronics. Currently, 
the majority of experimental measurements of conduc- 
tivity a in graphenei"— is analyzed on the basis the stan- 
dard Boltzmann approach within the Born approxima- 
tion predicting qualitative different results for short- and 
long-range impurity scattering^rJ^ 



a = 
a ~ 



const (short-range scattering), (la) 
n (long-range scattering), (lb) 



where n is the electron density. It has been recently ar- 
gued by Stauber et alr^ that a standard way to examine 
the collision rate within the Born approximation (utiliz- 
ing the unperturbed wave functions for a clean system) is 
not suitable for the case of short-range interaction such as 
vacancies, resonant impurities, cracks etc. Going beyond 
the Born approximation, Stauber et ali^ and Katsnelson 
and Novoselo\*i^ demonstrated that the short-range dis- 
order, with the accuracy up to logarithmic corrections, 
leads to a linear density dependency similar to the one 
for the long-range potential. 



o" = -; (ln-/7mi?o)^ , 

h Hi 



(2) 



where Rq is the scatterer's radius, and is the impurity 
concentration. (Similar results have been also obtained 
by Ostrovsky et ali^). Apparently, this has important 
consequence for interpretation of the experimental re- 
sults, as the linear density dependence of the conductivity 
is typically related to the long-range Coulomb impurities. 



and deviations from this dependence is attributed to the 
short-range scattering j^r— In contrast, Eq. implies 
the short- and long-range scattering may lead to similar 
density dependencies of the conductivity. Indeed, Mon- 
teverde et ali^ have recently analyzed the experiment 
on the basis of Eq. ^ and arrived to the conclusion 
that strong neutral defects (as opposed to the long-range 
Coulomb impurities) was the main scattering mechanism 
in graphene. The dominant role of neutral defects has 
been also recently outlined in Refs flTlflsl 

The reliability of the above predictions dH) , ^ can be 
established by testing them against "exact" Landauer- 
type quantum mechanical numerical calculations for the 
conductivity based on the tight-binding (or Dirac) Hamil- 
tonian for carriers in grapheneJ^"— . Recently, Adam et 
ali^ compared the standard Boltzmann and the Lan- 
dauer approaches for the case of a long-range Gaussian 
potential that varies smoothly on the scale of a lattice 
constant. 



Ni, 



= ^ Uii exp 



2e 



(3) 



where ^ can be interpreted as the effective screening 
length, and the potential heights is assumed to be uni- 
formly distributed in the range Ui € [—(5,(5]. The con- 
ductivity obtained in the two approaches agree quan- 
titatively away from the Dirac point, which was inter- 
preted as a proff of validity of both. According to Adam 
et ali^, the conductivity follows a density dependence 
a ^ r?!"^ . This however disagrees with all experimental 
observations reported so far— i^^"— and with previous 
Landauer-type numerical calculations22i2iiHi2^ demon- 
strating the linear or sublinear density dependence of the 
conductivity. Hence, a comparison between the Boltz- 
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mann and Landauer approaches still remains an open 
and an important issue. 

The main purpose of the present study is to com- 
pare the exact Landauer tight-binding (TB) conductiv- 
ities with those given by the standard Boltzmann ap- 
proach within the Born approximation as well as with 
those given by Eq. ([2]). As in the standard Boltzmann 
approach the density dependence is different for long- and 
short-range scatterers, one of our aims is to investigate 
whether the TB calculations also give different density 
dependencies for these scattering mechanisms. Finally, 
the Born approximation is valid for the case of weak scat- 
tering when the wave functions remains unperturbed. It 
is not however apparent the condition of a weak scatter- 
ing is satisfied in a parameter range typical for realistic 
systems. By comparing the exact TB calculations with 
those based on the predictions of the standard Boltzmann 
approach we test the applicability of the later to realistic 
systems. 
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FIG. 1: (a) The conductivity vs the electron density for two 
representative grapliene strips with different lengths L — 3000 
and 1500 sites for the case of the short-range scattering = 
a) with = 10 cm . Each curve is averaged over 8 

impurity configurations. The impurity strength 5 = 0.86t. 
(b) The inset shows the dependence 7 = 7(71). 



II. THEORY 

We calculate conductivity of the graphene sheets us- 
ing the standard p-orbital nearest-neighbor tight-binding 
Hamiltonian for non-interacting electrons for zero tem- 
perature, H = J2iVi N) (*l - N) O'L with the hop- 
ping integral t — 2.7 eVi^ We utilize a model for 
screened scattering centers of the Gaussian shape (Eq. 
([3])) commonly used in the literaturei^'^'^i^. The cor- 
relator of the potential (O has the form^i^^i^ (ViVj) = 

^2^-^^^ exp ^— ^'''^^^^^ ^, where the dimensionless im- 
purity strength is described by the parameter K « 
given by the screening length ^, 
the potential strength S, and the relative impurity con- 
centration Tiirnp (a being the carbon-carbon distance, and 
vp is the Fermi velocity). For realistic graphene samples 
K ^ 1 ^ 10, i9i20,24,25 ^ jj^^ pyj. calculations we use the 
screening length 1 < ^/a < 8, spanning the range be- 
tween the short-range potential = a), and the long- 
range potential = 8a) that varies smoothly on the 
scale of a lattice constant. Note that we also performed 
calculations for the (5-scattering, and the results obtained 
are, as expected, practically identical to those obtained 
for the case of the Gaussian disorder with ^ = a. 

The conductance G and the electron density n are 
computed with the aid of the recursive Green's function 
techniqu o^^i^*^'^^ . We assume that the semi-infinite leads 
are perfect graphene strips of the width W, and the de- 
vice region is a rectangular graphene strip W x L where 
the impurity potential is defined. The zero-temperature 
conductance G is given by the Landauer formula G ~ 
where T is the total transmission coefficient be- 
tween the leads. Then we calculated the conductiv- 
ity (7 = -^G, the electron density n = ^ dE D{E), 
the mobility /i = and the mean free path (mfp) 
mfp = as a functions of the Fermi energy Ep. 



The density of states (DOS) D{E) is computed by aver- 
aging the local density of states (LDOS) over the whole 
device area. The LDOS is given by the diagonal elements 
of the total Green's function.—. Because of computa- 
tional limitations we study the strips with L/W > 1 (In 
most calculations we use L = 368 nm (3000 sites) and 
L/W 6). We however checked that the obtained re- 
sults are insensitive to the aspect ration L/W as soon as 
L/W > 1. 

The results of the tight-binding calculations are com- 
pared to the predictions obtained within the standard 
Boltzmann approach within the Born approximation 

for the scattering potential (jS]). The conductivity of 

2 

graphene sheet is given by^"— a = e^rD^Ep)^, where 
the DOS D{Ep) = 2Ep/ (nfi^vj.) , the dispersion rela- 
tion E = hvpk, and k = y/nn, with vp = ^ being 
the Fermi velocity. The scattering rate within the Born 
approximation reads t-^ = ^D{E)1 /J' f 
where Uq is the Fourier transform of the scattering po- 
tential, q = |k — k'l = 2fcsin|, where 9 is the an- 
gle between the initial and final states k and k'. Us- 
ing the Wiener-Kitchine theorem for the correlator, we 
obtain the Fourier transform of the Gaussian potential, 
= A'(?iUi?)^ exp (— , which leads to the ex- 
pression for the conductivity^!^, 

_^ 4e2 7rn^2g7r<^ const, Trn^^ < I, , , 

where Ii is the modified Bessel function. 



III. RESULTS AND DISCUSSIONS 

The Boltzmann predictions for the density dependence 
of the conductivity are valid in the diffusive transport 
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FIG. 2: The conductivity of the graphene sheets vs electron 
density for the short-range Gaussian potential with £, = a for 
different potential strength S. Each curve is averaged over 
8 impurity configurations. Straight dashed lines show fitting 
a ^ rf . Red solid and dot-dashed lines show predictions based 
respectively on Eqs. Q and (b) The dependence 7 — 
7(^). (c) The m.f.p. vs the electron density n. 



FIG. 3: (a) The conductivity of the graphene sheets vs elec- 
tron density for the long-range Gaussian potential with ^ = 8a 
for different potential strength S. Each curve is averaged over 
8 impurity configurations. Straight lines show fitting a ^ rf . 
Dashed lines show predictions based on Eq. (j4|. (b) The 
dependence 7 = 7(5). (c) The m.f.p. vs the electron density 
n. 



regime when the mean free path is larger than a sys- 
tem size. Let us therefore first discuss a transition from 
the balhstic to diffusive regime focussing on the short- 
range scattering, ^ = a. In a purely ballistic regime (no 
impurity scattering) the conductivity follows the density 
dependence a ^ n"' with 7 = i.— It has been demon- 
strated for the case of the long-range Gaussian scatterers 
that with the increase of the system size the exponent 7 
gradually increases from its ballistic value reaching the 
value 7 = 1 in the diffusive regime.— Figure [T] shows 
the dependence a = a{n) for the case of the short-range 
scattering calculated within the TB approach. The con- 
ductivity shows the same behavior as for the case of the 
long-range scattering^^ with 7 increasing from ^ in the 
ballistic regime to 7 = 1 in the diffusive regime as the size 
of the system or the impurity concentration increases. 
This obtained density dependence in the diffusive regime 
(7 = 1) is in a stark contrast with the standard Boltz- 
mann predictions for the (5-impurity scattering, Eq. ()lap . 
when a is expected to be density independent (7 = 0). 

A more detailed comparison between the TB and 
Boltzmann calculations for the short-range scatterers for 
different impurity strengths S is presented in Fig. [5] (a)- 
(b). As expected, for very weak scattering {S < 0.5t) the 
transport is in the ballistic regime with 7 sa 1/2. This 
is fully consistent with the calculated m.f.p. which is 
comparable to the largest dimension of the system L, see 
Fig. [21(c). For the case of strong scattering {S > 2.5t) the 
system is in the diffusive transport regime when the cal- 
culated m.f.p. is smaller than the smallest dimension of 
the system W. In this regime the exponent 7 saturates 
to 1, see Fig. [2] (b). Figure [5] (a) also shows the con- 
ductivity calculated on the basis of the standard Boltz- 
mann approach, Eq. (U), as well as given by Eq. 
The Boltzmann theory predicts that 7 = which is in 
qualitative disagreement with the numerically calculated 



exponent 7 w 1. Boltzmann predictions arc also quanti- 
tatively different from the tight-binding calculations with 
^Boitz ^ ^TB (•j^Qi^g f^YyQ logarithmic scale of the figure). 
At the same time, we find that the TB calculations are 
in a good qualitative and even reasonable good quan- 
titative agreement with Eq. ^ predicting quasilinear 
density dependence of the conductivity. Why does the 
standard Boltzmann approach fail to describe the con- 
ductivity of the system at hand? Following Staubcr at 
alr^, we believe this is because the scattering rate r in 
the standard Boltzmann approach is calculated in the 
Born approximation, with unperturbed clean-graphene 
wave functions. Apparently, this approximation is appli- 
cable in the case of weak perturbations, but cannot be 
applied for strong scattering potential. In contrast, the 
approach proposed by Staubcr at al. uses wave functions 
for a hard- wall barrier, appropriate in the case of strong 
scattering. 

Let us now discuss the transition regime between the 
ballistic and diffusive behavior when W < m.f.p. < L 
(which corresponds to the impurity strength in the vicin- 
ity of (5 ~ t) where the exponent 7 shows a pronounced 
minimum dropping to 7 w 0.2 for 5 ~ t . We are not 
aware of any theories addressing this transition regime 
corresponding to the quasiballistic transport. We specu- 
late, however, that this peculiar behavior, with the con- 
ductivity becoming weakly dependent on the concentra- 
tion, might be related to the corresponding Boltzmann 
prediction of 7 = in the diffusive regime (even though 
the Boltzmann theory is not formally applicable in the 
case under consideration, with W < mfp < L). 

Let us finally compare the tight-binding and Boltz- 
mann calculations for the long-range Gaussian scatterers 
with ^ = 8a for different impurity strengths, see Fig. 
[31 As for the case of the short-range scatterers the TB 
calculations exhibit the ballistic behavior (m.f.p. ss L) 
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with 7 = 1/2 for weak scattering and the diffusive be- 
haviour (m.f.p. < W) with 7 « 1 for the strong scatter- 
ing. Again, the resuh obtained in the diffusive regime, 
7 « 1, is quahtatively different from the corresponding 
Bohzmann prediction. Indeed, for the latter case the 
exponent 7 is poorly defined because for the considered 
density interval tttt,^^ 1 which, according to Eq. (U), 
corresponds to the transition regime between two asymp- 
totes a =const and a ^ v?!"^. Besides, the Boltzmann 
and the tight-binding calculations disagree even quanti- 
tatively with 0"^°'*^ ^ (j^^ (Note the logarithmic scale 
of Fig. E] (a)). 

The opposite limit ttti^^ 3> 1 (when a ^ r?!"^) was con- 
sidered by Adam et al^ who found a good qualitative 
and quantitative agreement between the Landauer-type 
and the Boltzmann calculations. This regime (in contrast 
to the regime -Kn^ < 1 considered here) corresponds to 
high electron energies and smooth potential (with large 
^) when the scattering is weak and the Boltzmann the- 
ory within the Born approximation is therefore justified. 
However, the density dependence predicted by Eq. (|H), 
a ^ n^/^, has never been observed in any experiment. 
In contrast, the dependence a ^ n, obtained in the TB 
calculations for the regime im^^ < 1 considered here, is 
in agreement with the majority of experimental findings. 
We regard this as a strong indication that the regime ap- 
propriate for realistic graphene samples is Trn^^ ^ 1. As 
demonstrated here, in this parameter range the results 
obtained in the standard Born approximation disagree 
both quantitatively and qualitatively with those obtained 
by the exact TB calculations. This therefore questions 
the validity of the standard Boltzmann predictions within 



the Born approximation for realistic graphene sheets. 



IV. CONCLUSIONS 

(«) In relatively small systems the transport is in the 
ballistic regime for both short- and long-range scatterers, 
and the density dependence of the conductivity is cr '--^ 
712. As the system size (or the impurity concentration) 
increases the ballistic regime becomes superseded by the 
diffusive regime, in which the TB calculations predict the 
same linear density dependence, cr ~ n, for both short- 
and long-range scattering. 

(ii) In the case of short-range potential the obtained 
linear dependence is in quantitative and qualitative dis- 
agreement with the standard Boltzmann predictions 
within the Born approximation, Eqs. (|lap and (U), but 
in agreement with the predictions going beyond the Born 
approximation, 

(Hi) Even for the long-range Gaussian potential the 
standard Boltzmann predictions (0]) are in quantitative 
and qualitative disagreement with the TB results in 
the parameter range corresponding to realistic systems 
(^"-^^ ^ 1 regime). This questions the applicability of 
the predictions based on the standard Boltzmann theory 
for conductivity in graphene which are widely used for 
interpretation of experimental data. 
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